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Abstract 



I We present generic form of the covariant nonlocal action for infrared mod- 



ifications of Einstein theory recently suggested within the weak-field curvature 



J> ' expansion. In the lowest order it is determined by two nonlocal operators — 

■ kernels of Ricci tensor and scalar quadratic forms. In models with a low strong- 

■ coupling scale this action also incorporates the strongly coupled mode which 
, cannot be perturbatively integrated out in terms of the metric field. This mode 
! enters the action as a Lagrange multiplier for the constraint on metric vari- 
' ables which reduces to the curvature scalar and enforces the latter to vanish 

on shell. Generic structure of the action is demonstrated on the examples of 
p) ' the Fierz-Pauli and Dvali-Gabadadze-Porrati models and their extensions. The 

gauge-dependence status of the strong-coupling and VDVZ problems is briefly 
discussed along with the manifestly gauge-invariant formalism of handling the 
braneworld gravitational models at classical and quantum levels. 



1. Introduction: new mechanism of small cosmo- 
logical constant 

As is well known, the essence of the cosmological constant problem (CCP) consists in 
the enormous gap between the average density of energy in the modern Universe S ~ 
10~2^g/cm^ ~ (10~^eV)'^, generating according to Einstein equations the cosmological 
acceleration with the Hubble constant Hq, 

Hl^GE, (1.1) 

and the vacuum energy scale of all field theory models ranging from electroweak the- 
ory, £ ~ (ITeV)^, to quantum gravity, £ ~ (10^^ GeV)"^. Old attempts to resolve 
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this problem were based on biiilding the models with zero vacuum energy, incorpo- 
rating supersymmetry which protects the cosmological constant from renormalization. 
The cosmological constant mechanism in such models is based on cancellation of con- 
tributions of particles and their superpatners and stops working in the phase with 
spontaneously broken supersymmetry [1]. Moreover, it becomes completely meaning- 
less in the framework of the cosmological acceleration phenomenon [2, 3, 4] for which 
the magnitude of the cosmological constant, though being very small, is still different 
from zero. 

An alternative solution of CCP may be attributed to the curvature scalar sector 
of the Einstcin-Hilbcrt action, rather than to its cosmological constant term. The 
smallncss of in (1.1) can be explained not by small value of £^ but by the smallncss 
of the proportionality coefficient — gravitational coupling constant G. In other words, 
the cosmological acceleration is so small not because the vacuum energy is small, but 
rather because the latter is gravitating very weakly. Special property of the vacuum 
energy as compared to other local sources of the gravitational field is the degree of 
its spacetime homogeneity. It is assumed that it is not clustering and practically 
homogeneous at the horizon scale, 

^ ^0, (1-2) 

and gravitates with its long- wavelength gravitational constant Gld ^ Gp which is 
much smaller than the Planckian constant determining the everyday physics in the 
range of galaxies, planetary and solar systems, submillimeter Kavendish type experi- 
ments [5], etc.. Hi ~ Gld^ ^ GpS. 

This idea, that was apparently formulated in such explicit form for the ffist time in 
[6], represents the replacement of the fundamental gravitational constant in Einstein 
equations by the nonlocal operator, which for reasons of covariance can be regarded as 
a function of the d'Alembertian and which interpolates between the Planckian value 
of the gravitational constant and its long-distance magnitude^ 

G G{n), Gp > G{n) > Gld- (1.3) 

Note that this mechanism of scale dependent gravitational constant is not unique. 

The notion of scale includes not only the degree of spacetime inhomogeneity, but also 
the field amplitude. Therefore, the infrared modification of the theory can also be based 
on the gravitational "constant" locally depending on some distinguished physical fields 
— sort of quintessence [1, 8]. However, such modifications arc less universal because 

^The idea of scale dependent Newton constant was also put forward in [7], though it was not 
formulated in terms of the nonlocal operator. 
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they are attributed to the behavior of concrete quintessence field, while the mechanism 
of the nonlocal replacement (1.3) leads to modified Einstein equations 

- ^g^.uR = 87rG'(n) T^,, (1.4) 

in which independently of the field content of the matter source T^i, = T^i/(x) the grav- 
itational strength of the latter is determined solely by the degree of its inhomogeneity 
in X. 

One of the first problems associated with (1.4) is the lack of general covariance — 
no covariant action can generate such equations of motion with a nontrivial G{0). The 
solution of this problem was suggested in [9] by viewing (1.4) only as a first, linear in 
the curvature, approximation for correct equations of motion, their covariant action 
being constructed as a weak field expansion in powers of the curvature with nonlo- 
cal coefficients. The most solid example of the theory incorporating the mechanism 
(1.3) is the Dvali-Gabadaze-Porrati (DGP) model of brane-induced gravity [10] which 
is very interesting due to the fact that it naturally contains the mechanism of the cos- 
mological acceleration [11]. However, this model turned out to be suffering from the 
strong-coupling problem [12, 13, 14] which invalidates the curvature expansion theory, 
makes it inefficient without the fundamental UV completion at the quantum level and 
features the van Dam-Veltman-Zakharov (VDVZ) discontinuity problem [15]. Though 
this problem does not indicate the physical inconsistency of the underlying theory and 
is likely to be circumvented in tree-level applications by accounting for nonlinear effects 
[16] or using the relativistic expansion [17], in quantum domain it rises at full height, 
because of the essentially perturbative nature of Feynman loop expansion. Therefore, 
its status becomes very important for possible generalizations of the Einstein theory. 

The goal of this paper is to discuss the most general version of the covariant nonlo- 
cal action for long-distance modifications of Einstein theory, suggested in [9] . It turns 
out that in the lowest (quadratic in curvature) order it is determined by two nonlocal 
operators, kernels of Ricci and scalar curvature quadratic forms, and also may contain 
additional term responsible for the strong-coupling problem. The strongly coupled 
mode enters this term linearly like the Lagrangian multiplier of the auxiliary con- 
straint on metric field and cannot be perturbatively integrated out in terms of metric 
variables. After discussing this in Sect. 2 we demonstrate generic structure of the action 
on the examples of Fierz-Pauli and DGP models, consider the details of VDVZ and 
strong-coupling problems and review their generalizations which allow one to circum- 
vent these problems. In particular, in Sect. 5 we review a recently suggested constrained 
perturbation theory in the DGP model [18] and show that it turns out to be a nonco- 
variant modification of the original DGP model. In the concluding section we briefiy 
discuss the gauge independent status of the VDVZ and strong-couphng problems and 
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outline the manifestly gauge-invariant technique for classical and quantum braneworld 
models. 



2. Covariant nonlocal action for infrared modifica- 
tions of Einstein theory 

The idea of replacing the gravitational constant by a function of the covariant d'Alember- 
tian □ = Sf^^VoV/j, according to [6] consists in the following modification of the left 
hand side of Einstein equations 

R^u- 7;9t.uR] = 7;Tf,u, (2.1) 



167r V 2"^'' 7 2 

where the nonlocal Planck mass is a function of the dimensionless combination of □ 
and additional length scale L, interpolating between the Planck constant for matter 
sources of small size <^ L and long distance constant Gld = G{0), 



G{n) 

If the function of z = L'^O satisfies the conditions Tiz) ^ at 2; ^ 1 and J-'{z) — >■ 
JF(0) ^ 1 at 2; ^ 0, then the infrared modification is inessential for processes varying in 
spacctimc faster than 1/L and vicy versa — becomes significant — for slow phenomena 
with wavelengths ~ L and higher. 

An obvious difficulty with the construction of the above type is that for any nontriv- 
ial operator T{L'^U ) the left hand side of (2.1) does not satisfy the Bianchi identities 
and cannot be obtained by varying the covariant action. In particular, a naive attempt 
to modify the gravitational action according to 

Ml jdxg^/^R^ J dxg^/^M\n)R = M\0) jdxg^^^R (2.3) 

is meaningless, because after the integration by parts the covariant d'Alembertian when 
acting to the left (on 1) picks up its zero mode and the nonlocal operator in all regimes 
reduces to its infrared value Afp(O). 

This problem can be solved within the assumption of weak field approximation 
implying that the equation (2.1) represents only the first linear term of the perturbation 
expansion in powers of the curvature. Its left hand side must include higher orders 
in the curvature, the nonlocal gravitational action <S'jvl[5';ui/] generating the modified 
equations according to 



5SNL[g] ^ M^u) 



l^..^^+o[i?JJ. (2.4) 
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To obtain the leading term of SNLldfiu], this equation can be functionaUy integrated 
in the cxpUcit form [9] with the aid of the covariant curvature expansion technique of 
[19]. The essence of this technique consists in the possibihty to convert noncovariant 
series in powers of gravitational perturbations h^i, into the series of spacetime curvature 
and its derivatives with the covariant nonlocal coefficients. The starting point is the 
expansion of the Ricci tensor 

i?^, = -1 □ V + ^ (V^^. + V.F^) +0[hl,], (2.5) 

{Fn = V^/t^A ~ I Vj^/i is the linearized de Bonder- Fock or harmonic gauge), which can 
be solved by iterations with respect to h^^y in the form of the nonlocal expansion in the 
curvature, beginning with 

V = -5^M- + ^^^^ + V./m + 0[R%]. (2.6) 

Here V^fi, + Vufn reflects the gauge arbitrariness in this solution originating from the 
terms with a harmonic gauge in (2.5). 

As a result, the nonlocal action generating (2.4) begins with the quadratic order in 
the curvature [9] 

SnlM--^ J dxg'/'l^(R^-'-^g^''R)^^^R,^ + 0[R%]y (2.7) 

Interestingly, in the simplest case of constant M^(n) = Mp it should reproduce the 
Einstein- Hilbert action, which looks puzzling because it does not at all contain the term 
linear in the curvature. The explanation of this paradox consists in the observation 
that the Einstein action in the asymptotically-flat spacetime with the metric behaving 
^ gnu — ^nv + ^/ii/, — O (l/|xp), \x\ oo, includes the Gibbons-Hawking surface 
integral over spacetime infinity 

SE[g,A = -Tf dxg'I^R{g) + -f / da^ {d'^h,, - d,h). (2.8) 

This integral can be transformed to the form of the volume integral of d^{d''h^i, — d^h) 
— linear in /i^^ part of the scalar curvature — and, with the aid of (2.6), covariantly 
expanded in powers of the curvature. Up to quadratic terms inclusive this expansion 
has the form [9] 

j^^^^Ja^id^K, - d,h) = j dxg'/' {i? - (r^'' - + ...| . (2.9) 

Therefore, on substituting to (2.8) linear in Ricci scalar terms cancel out, and the 
quadratic terms reproduce the expression (2.7) with the coefficient M^(n) = Mp which 
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can be pulled out of the integral. The result is the nonlocal form of the local Einstein 
action [20, 21, 9], 

SE[g,.] = -^J dxg'/'l^{R'^-^ - \g^^R)^R,^ + 0[4j| . (2.10) 

The fact that this action begins with the quadratic order in R^j^^, ~ hj^^ obviously 
corresponds to the theory of massless spin-2 field. Less trivial is the nonlocahty of this 
action, which is the price one should pay for the manifest covariance of this expansion 
in contrast to the local in terms of hn„, but noncovariant action for symmetric spin-2 
tensor field. 

Thus the action (2.7) serves as a direct realization of the idea of nonlocal gravita- 
tional coupling constant — replacement of Mp in the nonlocal version of the Einstein 
action (2.10) by the operator (2.2) as suggested in [9]. However, there is a question 
(that was not addressed in [9]) to what an extent this infrared modification is generic 
even in the quadratic order approximation? It is obvious that in the general case the 
operator kernels in the quadratic Ricci tensor and scalar forms can be different, so that 
the generalization of (2.7) takes the form 

SML[g,A --^J dxg'/' -lR^^R + 0[Rl]j , (2.11) 

where the two nonlocal "masses" should tend to one and the same Planckian limit Mp 
only in the short and intermediate distance domain, 

Ml^ia) Ml, □ > 1/L. (2.12) 

The last requirement implies the recovery of the Einstein general relativity at in- 
termediate distances including, in particular, the absence of the VDVZ discontinuity 
problem — correct tensor structure of the gravitational potential. For the theory (2.11) 
the linear potential generated by the conserved matter sources T^^, has the form 

V = - .ifl^ (t,. - ^ «(□) t) + d,i. + d.C,. (2.13) 



M2(n)n V 2 

Here the longitudinal terms correspond to gauge arbitrariness and the operator coeffi- 
cient «(□) equals 

2M|(D) - (□) 

= 3M|(n)-2Mf(n)- ^'-'^^ 

Generically, it should take the general rclativistic value a = 1 only in the limit 
(□) ^ M|(n) accounting for the restriction (2.12). 
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Expression (2.11) seems to give the most generic form of the gravitational theory 
in weak field approximation, dictated by general covariance. However, this is also not 
the end of the story — infrared limit of the theory together with spin-2 polarization 
can accommodate additional degrees of freedom that are not taken into account in this 
expression. One would think that, under a natural assumption that these degrees of 
freedom are not directly coupled to matter, they can be integrated out, which reduces 
to additional contributions to M^{0) and M|(n). However, such a reduction is always 
possible except the case when the equations of motion for additional fields cannot be 
solved for them in terms of the metric. This happens when these fields enter the action 
linearly and play the role of Lagrange multipliers of certain combinations of metric 
variables. Unfortunately, such fields persist at higher orders of perturbation theory 
and, thus, give rise to the problem of the low strong-coupling scale. This precludes 
from using the conventional perturbation theory in a wide range of distances very 
well probed by numerous Kavendish type and solar system experiments. Below we 
demonstrate this phenomenon on the examples of the Fierz-Pauli (FP) theory and the 
effective 4-dimensional theory of brane induced gravity of Dvali-Gabadadze-Porrati 
(DGP). 



3. Fierz-Pauli model and VDVZ problem 

The simplest infrared modification of the Einstein theory is represented by the Fierz- 
Pauli model of massive tensor field. It is described by the quadratic part of the Einstein 
action (2.8) modified by the noncovariant mass term on the /Zat-spacetime background 

^.ass[^..] = / ( X - T ' 

V = 9^1^ - Vfiu: h = rj^'^hij,^. (3.2) 

This is the only Lorentz-covariant combination of mass terms which guarantees the 
absence of ghosts in the theory. Linear equations of motion in this model have the 
form 

1 Tfl^ 

Rixu - ^Vi^i'R + ^ {h^^ - f]i,^h) = SttG^ Ti^^, (3.3) 

where under i?^,^ we understand the part of Ricci tensor linear in the graviton field. 
Differentiating this equation and taking into account the linearized Bianchi identity 
and the concservation of we obtain the "gauge" for h^i, 

9''(V-^/.-^) = 0, (3.4) 
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which as a corollary has the vanishing of the linearized curvature scalar 



R = d^d'^hf,^ -nh = (3.5) 

(note that this equation is satisfied even for nonvanishing trace of matter stress tensor). 
As a result the gravitational field generated by the matter source takes the form 

V = -167rG4 (t^. - lvf..T^ + d^Cu + d,^^ (3.6) 

up to the longitudinal terms which do not couple to the conserved matter sources. Since 
the FP theory is not gauge invariant, these terms are not arbitrary and determined by 
the fixed vector 

«-.= -?^^9,r, (3.7) 

which, in particular, guarantees the validity of Eq. (3.4). 

Tensor structure of massive graviton field here differs from the general relativistic 
case — the trace of T^,y in Einstein theory has a coefficient 1/2 rather than 1/3 in FP 
model. This discrepancy remains also in the limit of vanishing mass and constitutes the 
VDVZ discontinuity problem [15]: massless hmit of FP model does not reproduce the 
predictions for massless graviton of Einstein theory. This problem originates from the 
additional degree of freedom which is missing in Einstein theory and, on the contrary, 
exists in FP model for all values of the graviton mass. From the viewpoint of the general 
framework for infrared modifications of Einstein theory, this degree of freedom should 
enter the action (2.11) as a Lagrange multiplyer responsible for the auxiliary equation 

(3.5) . Without it the equations of motion following from (2.11) do not reproduce the 
tensor structure (3.6) with any choice of operators (□) and M|(n). Indeed, the 
coefficient (2.14) takes a numerical value of FP model, app — 1/3, only in a singular 
limit (□) — > 0. On the contrary, the inclusion of the constraint (3.5) in the action 
(2.11) with a Lagrangian multiplier — additional scalar field — improves the situation. 
It leads to the equations of motion which yield as a solution the gravitational potential 

(3.6) for the following choice of nonlocal operators 

1 □ — 

Ml{U) = Ml{U) = —^^. (3.8) 



4. DGP model and the strong coupling problem 

Effective long-distance modification of the Einstein theory is also provided by the DGP 
model of brane induced gravity. Its fundamental action includes the five-dimensional 
(bulk) Einstein term of the metric Gab{X), A = 0, 1, 2, 3, 5, and the Einstein term of 
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the metric g^y{x), // = 0, 1, 2, 3, on the 4-dimcnsional brane viewed either as a boundary 
of the 5D bulk or as a fixed point of Z2-identification on the 5D orbifold 

Sdgp 

[ Gab {X) ] = -J— [ d^X G^''^ R^^^ ( Gab ) 

The brane term is accompanied by the Gibbons-Hawking surface integral of the trace of 
extrinsic curvature of the brane [K] (analogous to the one at infinitely remote boundary 
written in (2.8) in the noncovariant form). The 4-dimensional gravitational constant 
here G4 — 1/M| is essentially different from the 5-dimensional one G5 = They 
both determine two different scales of the DGP model — M4 typically of the Planckian 
value and the crossover scale to the infrared regime 

(usually identified with the Hubble scale of the Universe m — H — 10~^^ cm~^). 

As is known, the DGP model is analogous to the FP theory with the mass term 
more soft in the infrared, m^/—D, in which the role of m is played by the scale 

(4.2). To show this, one can build the effective brane action which follows from (4.1) by 
integrating out the bulk metric Gab{X) subject to fixed boundary data on the brane. 
This can be done perturbatively. For this purpose expand the fundamental action in 
gravitational perturbations 

GABiX)^riAB + HABiX), (4.3) 

solve the linear equations for Hab{X) in the bulk and substitute the result back into 
the quadratic part of the action. We shall work in the coordinate system in which the 
brane is located at the fixed value of the fifth coordinate, = y — 0. For gauge fixing 
in the bulk add to the action (4.1) the gauge breaking term quadratic in the linearized 
de Donder-Fock gauge 

W[ ^ab] = - 2 y "^'^ V^^'FaFb, (4.4) 
Fa^O^Hab-^OaH. (4.5) 

In this gauge the bulk equations of motion are most simple and form the following 
boundary value problem 

n,HAB{X) = 0, (4.6) 
Hab{x, y) \y^Q = hAB{x), hAB{x) = ( ^n(^)^ hsix)). (4.7) 



9 



On a flat background Ds = □ + 9^, and the solution of this problem nonsingular at the 
bulk infinity can be written down in the following elegant form [13] 

HAB(x,y)^e-y'^hAB{x) (4.8) 

in terms of the auxiliary operator 



A = (4.9) 

(the case of Lorentzian spacetime we treat as the analytic continuation from the Eu- 
clidean space where the operator □ is negative definite, so that HAB{x,y) vanishes at 
y oo). 

Substituting the obtained solution in the 5D part of the DGP action (4.1) S^[Gab] 
(with the fivc-dimcnsional curvature and the Gibbons-Hawking integral) and taking 
into account the gauge breaking term (4.4) one gets [13] 

S^[Gab] + W[^ab] = -[^J J (^-/I'^'^AV + ^hAh 

+ hAh55 - ^/issA/iss) , (4.10) 

where m — 2M^/M^ is the DGP scale (4.2) and is the following combination of 
the metric induced on the brane and shift functions in the fifth dimension G^^ — 

V = V + ^ i9>.N, + 9uN^)- (4-11) 

Note that this combination is gauge invariant with respect to the four-dimensional 
coordinate transformations 

S^K. = d^C, + d,C^, 5^N^ = -A^^, 5^h,5 = 0. (4.12) 

In their turn, these transformations represent the restriction to the brane of the residual 
gauge transformations in the bulk 6~Hab = c^^Ss -|- Ob'^a, 

E^^ix, y) = e-^^ e(a;), El'{x, y) = 0, (4.13) 

which, because of □5S^(x, y) = 0, leave invariant the de Bonder gauge (4.5) and do 
not shift the brane from its location at y = 0, . 

Thus, as expected, the effective action (4.10) induced from the bulk is invariant with 
respect to four-dimensional transformations, but this invariance is actually realized by 
means of the Stueckelberg fields A*"^, usually introduced by hands for the covariantiza- 
tion of the noninvariant action. By varying (4.10) with respect to A'^, these fields can 
be excluded in terms of the metric variables 

= i V - \9,h - \d,h,^ , (4.14) 
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which leads to the manifestly invariant expression for h/j^i, in terms of the linearized 
Ricci tensor 

1 d d 

V = -'2-R,u + -^hs. (4.15) 

Their substitution to the bulk action leads to 

^5 + -^gauge = ^ ^ jd'^x (-R'^^^R^, + - R^h^) , (4.16) 

where the scalar variable h^^ can no longer be excluded in terms of the metric field — 
the situation discussed above. 

Adding (4.16) to the four-dimensional part of the DGP action, rewritten in the 
nonlocal form (2.7), we finally obtain the quadratic part of the effective action on the 
brane 

^LGp[^..,n] = dx^V2|(^.._ i^..^)£_^^^^ + ^ni?|, (4.17) 

where the Lagrangian multiplier to the scalar curvature 11 is related to h^^ component 
of the 5D metric according to 

/i55 = -2An. (4.18) 

The variable 11 was introduced in [13] as a longitudinal part of the 5D shift function 
— di^Il + N'^. It parameterizes the brane bending in the form of the 5D diffeomor- 
phism of the bulk metric 

6^H,, = d^E^, (4.19) 
5si/55 = 2 dyE^ = -2An, (4.20) 

with the vector field 'EA{x,y) — 5^ e~^'^n(x). Similarly to (4.13) this diffeomorphism 
does not break the de Bonder gauge in the bulk, but it shifts the brane in the direction 
of the fifth coordinate by S^(x, 0) = n(x) and, therefore, is not the symmetry of the 
action. It does not manifest itself in the bulk, and its effect reduces to the contribution 
on the brane, which begins with the local term mlli?. 

Thus, the 4-dimensional effective action of the DGP model does not take the form 
(2.11), because it contains auxiliary constraint with the Lagrange multiplier which 
cannot be expressed in terms of metric variables. Therefore, in the linear approximation 
the DGP model on the brane is effectively described by the Fierz-Pauh theory with 
the nonlocal mass term of the form (4.16), generated from the bulk. In essence, the 
expression (4.16) turns out to be the covariant completion of this term (3.1) with the 
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infrared soft mass ^ V niA. Linearized gravitational potential of matter source in this 
model is analogous to (2.13) 




Interaction of matter on the brane is determined here by the propagator which for 
small spacetime intervals |x| <^ L = 1/m, |n| 3> m^, obviously coincides with the 
4-dimensional one. On the contrary, beyond the crossover scale (4.2) the coupling 
becomes five-dimensional. This is usually interpreted as a gravitational leakage to the 
bulk — 4D graviton is metastable and decays within a lifetime L = 1/m similarly to 
the Gregory- Rubakov-Sibiryakov model [24]. In contrast to GRS model which suffers 
from ghost negative-energy states [25, 26, 27], the DGP model (hke FP theory) is 
ghost- free. 

Similarly to FP model, the constraint term of (4.17) results in the VDVZ problem — 
tensor structure of the gravitational potential (4.21) differs from general relativity and 
corresponds to to Ficrz-Pauli theory in all ranges of distances. Note that kincmatically 
the variable 11 of brane embedding into the bulk is analogous here to the radion mode 
in the Randall- Sundrum model [28], which guarantees there the recovery of the general 
relativistic structure. Here, however, this mode does not cope with this task and the 
DGP model suffers from the VDVZ problem [10]. 

Another consequence of the constraint term in the DGP action (4.17) is the problem 
of low strong-coupling scale. Point is that the variable 11 has a nature of the Lagrange 
multiplier only in the quadratic order, while in higher order terms of the weak field 
expansion it enters the action nonlinearly and gives rise to composite operators built 
of the powers of 11, N^j,, h^^ and their derivatives. On the other hand, its kinetic 
term in the quadratic order originates entirely from mixing with h^^iy (in view of its 
Lagrangian multiplier nature in MpmUR ~ MpmdUdh) and is small because of 
smallness of m. As a result, after the diagonahzation of the full quadratic term the 
11- mode acquires the kinetic term ~ M|,m^ (511)^, and the transition to the canonically 
normalized field II, 11 = II/ (mMp), when expanding in 11 gives rise to higher and 
higher negative powers of small quantity m [12]. Then, the composite operators of 
high dimensionalities become suppressed by the factors of the form 1/Mpm'^ and get 
strong at the scale Ap^g = {Mpm'^y^^'^'^K As shown in [13, 29], the lowest scale occurs 
for the cubic in 11 interaction 



^Covariant structures of such a type as a realization of the nonlocal cosmological "constant" were 
also discussed in context of the renormalization theory in [22, 23]. 




(4.22) 
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and its value 

A = (m^Mp)^/^ ~ (1000 km) (4.23) 

is much lower than the submillimctcr scale (0.02mm)~^ up to which the Newton law is 
verified by high precision table-top experiments [5] . Analogous situation is well known 
for the nonlinear FP model (with the full Einstein term and quadratic mass term) — 
its strong-couphng scale being equal to A5 = (m^Mp)-*^/^ [12]^- 

The recovery of Einstein phase in the DGP theory (with the correct tensor structure 
of the propagator) can be attained by its synthesis with the Randall-Sundrum model. 
If the brane has a positive tension fine tuned to the negative cosmological constant in 
the bulk, then in the perturbative domain below the strong coupling scale, □ <C A^, 
the gravitational potential (2.13) has operator functions [31] 



M^{n) = Ml 
a{n) = 



1 + 



AK2(/A) 
2 Im 



2 + ml 3 2 + lm 



1 



(4.24) 



(4.25) 



IAK2{IA) 

in terms of McDonald functions of the first and second order, iTi 2(0;), and curvature 
radius / of the AdS bulk. In the distance range 1/m ^ 1/A ^ / this potential 
describes the 4-dimensional general relativistic law with «(□) — 1 and effective Planck 
mass Mp = M| (1 + lm/2) ~ M| [31]. Thus, within the hierarchy of the horizon 
(DGP) and the AdS scales 1/m = 1/if > Z this model does not suffer from the VDVZ 
problem. Actually, this is a generahzation of the well-known result that this problem 
is absent for spin-2 massive field in (A)dS spacetime [32, 33] with the cosmological 
constant A in the limit m^/A — > 0. 

Another generalization that helps to resolve the VDVZ problem consists in increas- 
ing the number of extra dimensions in the DGP model, which has a good motivation 
from string theory [34]. Correct tensor structure in this case is, however, achieved by 
the price of ghost states of tachyon nature [26] which indicates classical and quantum 
instabilities. Two approaches were suggested to avoid these states, based on the ne- 
cessity to make the UV regularization for branes with codimension N > 1 [35]. The 
regularization smearing the brane was used in [36] (together with the D — A + N- 
dimensional scalar curvature replacing on the brane the 4-dimensional one) and the 
model was shown to be ghost and tachyon free, and, moreover, its 5D version acquiring 
a big strong-couphng scale Ag = (Mjm^)^/^ :$> A3 much exceeding (4.23) [37]. 

■^At the classical level the strong-coupling and VDVZ problems in the FP model may apparently 
be circumvented by infinite resummation of nonlinear terms [16]. However, at the quantum level the 
theory has the strong coupling scale [30] which at best can be raised to (m^Mp)^/'^ by the inclusion 
of higher-dimensional operators [12]. 
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Another approach to DGP models with N >2 [38] was based on the interpretation 
of Green's functions poles different from [26]. The unitarity of the theory was recovered 
by noting that the ghost tachyon of [26] belongs to the unphysical sheet of complex 
Mandelstam variable and, apparently, can be removed by the sort of the Lee- Wick 
prescription in local field theory with ghost states [39]). The accompanying loss of 
analyticity of the propagator results in the loss of causality at the horizon scale L — 
Mi/M'^, which as advocated in [38] does not contradict causality at small and interme- 
diate distance scales. 

5. Constrained DGP model 

VDVZ and strong coupling problems in the DGP model invalidate the weak field ex- 
pansion everywhere except the ultrainfrarcd domain of energies E below the horizon 
scale E <^ m = 1/Hq. This makes the fiat-space perturbation theory (and its covari- 
antization according to the lines of Sect. 2) inefficient. Moreover, it turns out to be 
incomplete without the the UV completion at the quantum level, because the infinite 
tower of higher dimensional operators limits its predictive power even for larger dis- 
tances. Recently there were several attempts to save the situation [29, 18]. One of 
them [29] is based on the transition to the dynamically more safe non-flat background 
on which the nonvanishing brane extrinsic curvature K^^, introduces big kinetic term 
for the brane bending mode and, therefore, acts like a dilaton field controlling the 
strength of its interaction at each spacetime point. More radical is the suggestion of 
[18] to interpret the strong coupling problem as a gauge artifact that can be circum- 
vented by a special modification of perturbation theory called in [18] the constrained 
perturbation theory. 

This modification consists in adding to the action (4.1) the full gauge- fixing term 

which includes the set of special 5-dimensional gauges in the bulk 

B^ = d,H,,-d^H^f„ B^^d'^H^,, (5.2) 
along with the 4-dimensional harmonic gauge for the intrinsic metric h^j^^, on the brane 

i^M = 5' V - l^^h. (5.3) 

Prom the viewpoint of gauge-fixing procedure the brane gauge term is redundant, 
because the bulk term in (5.1) alone completely fixes all coordinate invariances of 
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the theory^. Therefore this term was interpreted in [18] as an additional constraint 
extending beyond conventional gauge-fixing procedure. The calculations show that 
for particular limiting values of " gauge- fixing" parameters ^ 0, 7 ^ and a = 1 
the gravitational potential of matter sources on the brane does not suffer from strong 
coupling and VDVZ problems. In particular, in Eq.(3.6) it has the operator coefficient 
«(□) interpolating between the FP value app — 2/3 in the deep IR regime A = 
V— C <^ m and the Einstein value a = 1 for short distances A ^ m and does not have 
terms singular for m — > 0. 

This result, however, can hardly be interpreted as a proper solution of these prob- 
lems in the DGP model. Rather than being the constrained perturbation theory for 
this model, the suggestion of [18] represents the modification of the model itself, and 
this modification is not unique. Indeed, the theory with the action (4.1) modified by 
the term (5.1) in the limit 

/?^0, 7^0, (5.4) 

but with an arbitrary parameter a, has the following gravitational potential on the 
brane, Hab = Hab\, (in what follows the vertical bar denotes the restriction of a bulk 
quantity to the brane) 

V = -M (^"^ - i-C),.r) . ^a.a„^T, (5.5) 

_ 167r 1 n-2amA 1 

M| 2n-4mA-3am2 □ ' ^ ' 

h^. = (5.7) 

with the operator coefficients (see the derivation in Appendix) 

□ - 3mA - 2am^ , , 

aa(n = ^, (5.8) 

' □-4mA-3am2' ^ ' 

The data (5.5)-(5.7) is propagated from the brane into the bulk by the massless Klein- 
Gordon equation^ (4.6) according to (4.8) and does not have dangerous parts blowing 
at m — > 0. 



^The residual gauge transformations for the bulk gauge (5.2), 6Hab = 29(^2^), are locally 
parametrized by one scalar field S5 propagating in the bulk, DsSs = 0, = 9;^(l/n)95S5, but 
the corresponding diffeomorphism should not move the brane, S5I =0. Therefore, with this condition 
on the brane and zero boundary conditions at infinity S5(X) vanishes identicallty everywhere in the 
bulk and thus does not leave any gauge freedom. 

^This is a corollary of the limit (5.4). For nonvanishing f3 and 7 only the transvcrse-tracclcss, -ffJJ", 
and conformal, r/^i/^, sectors of Hab are subject to this propagation law (see Eqs.(A.13)-(A.15) in 
Appendix). 
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The expression (5.5) for on-brane metric can be obtained from the effective action 
of the form (2.11) with the nonlocal form factors 

This action, in contrast to (4.17), does not contain any constraint term, because no 
strong coupling mode is present and all variables express in terms of the 4-dimensional 
metric. Since the constrained theory is not gauge- invariant, the longitudinal term of the 
solution (5.5) is not ambiguous just like in the Fierz-Pauli model with (3.7). Therefore, 
this solution can be regarded as obtained from the covariant effective action (2.11) in 
a particular (j-dependent gauge 

a-V-l °-'^^> ,/.^o. (5.11) 

This gauge interpolates between the Fierz-Pauli gauge (3.4) in the infrared and the 
harmonic gauge for high energies, just like the operator coefficient (5.8) interpolates 
between the corresponding Fierz-Pauli and Einstein values. 

Thus, for generic nonsingular a the brane metric (5.5) has the Einstein high-energy 
limit and the (soft-mass) Fierz-Pauli behavior at low energies. But at intermediate en- 
ergies it corresponds to the one-parameter family of physically inequivalent theories. In 
particular, they are not equivalent to the original DGP model, because their equations 
are spoiled by gauge terms which arc nonvanishing even on shell. For example, one 
of the bulk Einstein equations in the DGP model implies the vanishing 4-dimensional 
scalar curvature, ~ i?*^^) = 0, while for (5.5) it reads 

/?(4)__^ \ rr (r.^2\ 

M| □ - 4mA - 3am2 ^''■^^> 

and goes to zero only in the singular limit a ^ oo (this limit obviously corresponds to 
the original DGP model and is marred by the strong-coupling and VDVZ problems). 

The constrained theories depend not only on the cr, but also on the bulk "gauge- 
fixing" parameters (5 and 7. This gauge dependence is an artifact of introducing the 
redundant brane "gauge" in (5.1). As shown in Appendix, the Israel junction condition 
which underhes the derivation of (5.5)-(5.7) estabhshes the relation (A. 23) between this 
gauge and the boundary value of the bulk gauge \ . Only their combination (A. 23) 
is vanishing on shell, so that separately Ffj_{x) and Ba{X) are nonzero^. Therefore these 
constraints go beyond a consistent gauge-fixing procedure and modify the original 
model. 



^Brane gauge -F^ can be separately put to zero by demanding the continuity of -ffs^-coefficients 
across the brane [40], but this requirement seems contrived, because it does not follow from the geo- 
metrically invariant formulation of the Z2-symmetry of brane orbifolds. 
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6. Discussion and conclusions 



The above discussion shows that the strong-couphng and VDVZ problems are very 
robust within covariant infrared modifications of Einstein theory. In particular, for 
DGP model the VDVZ problem is a direct corollary of the 5-dimensional Einstein 
equations. The (55)-component of these equations implies the vanishing of the 4- 
dimensional scalar curvature (3.5) which means that there is no conformal term ri^j^^ip in 
the transverse-traceless decomposition of the metric potential on the brane h^iy = hj^^ + 
d^^v + du^^. Therefore, up to gauge transformation terms the linearized gravitational 
field is given by the transvcrsc-tracclcss part of the matter stress tensor, h^^ ~ T'^J i 
and this unambiguously fixes the a-coefficient of (2.13) at the Fierz-Pauli value 2/3, 
because 

TlJ = T,. - \r),.T + \d,d^^ T. (6.1) 

Thus, the VDVZ problem in DGP model cannot be circumvented without spoihng the 
basic dynamical equations in the bulk. This is exactly what happens in the constrained 
DGP model of [18]. 

The strong coupling problem also has a gauge-invariant origin, even though the 
derivation of Sect. 4 seems to point at the particular harmonic gauge as a source of this 
difficulty. Note that in the derivation of (4.17) we could not exclude via the variational 
equation the h^^ metric component which turned out to be the strongly coupled mode 
(4.18) — the Lagrangian multiplyer of the scalar curvature constraint in (4.17). The 
impossibility to solve the linearized Einstein equations in the bulk for lapse and shift 
functions Na — G^a follows from the degeneracy of the operator 

On flat space background it is really degenerate because its (5A)-column is identically 
zero and the (/xi/)-block, S^" — d^d'^ — ri^^O^ has as a zero mode an arbitrary longitu- 
dinal vector — di,Il — the (5i/)-component of the strong couphng mode (4.19). 

At this point it is useful to compare the gauge (in) dependence properties of the 
original DGP model and its constrained version of Sect. 5. The dc Dondcr gauge Fa, 
(4.5), used in Sect. 4 is qualitatively different from the gauge Ba in (5.2). It belongs 
to the class of relativistic gauges for which the residual transformations are locally 
parametrized by the full set of five propagating modes, DsS^ = 0, and the Faddeev- 
Popov operator Ds is dynamical for all components of S^. The fifth component S5 is 
ruled out by Dirichlet boundary conditions on the brane S5I = (because the movement 
of the brane is not a gauge transformation), and the residual transformations reduce 
to 4-dimensional diffeomorphisms (4.13). Their fixation requires four extra gauges 
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imposed on the branc. and this is a main point of departure from the special gauge 
(5.2) of [18]. In contrast to [18] these extra gauges do not violate Ward identities and 
do not change the physical contents of the model. This property is based on the fact 
that, though the effective brane action (4.10) is gauge dependent, its gauge dependence 
gets eradicated by the transition to the effective action (4.17) with on-shell values of 
the functions Na- This is different from the constrained version of the DGP model 
for which the effective form factors (5.10) explicitly depend on all - bulk and brane - 
"gauge-fixing" terms. 

To clarify the properties of the action (4.10) versus (4.17), note that the former was 
obtained as a fundamental DGP action 

S'^Gp[9,^, Na] = SogpIGab] L,,=G..[.,.,iv.] (^-^^ 
g^.u = G^,\, Na = G5a\, (6.4) 

evaluated on the solution of the bulk Einstein equations GAB[9fiu, Na] subject to inde- 
pendent boundary conditions for all five-dimensional metric coefficients on the brane. 
The action functional for such boundary conditions depends on the choice of gauge, 
because its on-shell restriction is not complete regarding the boundary values Na{x) 
of G5^(X)-coefficients. Only their exclusion by virtue of variational equations 

in terms of the metric g^,^, Na — NaIquu], makes the action gauge independent^. 
However, the strong-coupling problem prevents from solving these equations for all 
A^^ = NaIo^u] ^ only its transverse yU-components can be excluded, and at least one 
scalar mode 11 remains off shell among the arguments of (4.17). Nevertheless, the 
result turns out to be covariant and bulk-gauge independent, though suffering from 
the strong coupling and VDVZ problems. 

In essence, the gauge difficulties with these problems follow from the lack of mani- 
festly covariant and gauge independent formalism. The on-shell reduction which makes 
the brane action gauge independent does not resolve this problem. In particular, the 
construction of brane-to-brane propagator requires the off-shell extension of the brane 
action and, therefore, is sensitive to the choice of gauge. Moreover, the off-shell ex- 
tension of the action is necessary at the quantum level, where it should be integrated 

^In the rclativistic gauge Fa{X), (4.5), the variational equations (6.5) represent the requirement 
of zero boundary data for these gauge conditions, Fa\ =0. In the bulk they satisfy the homogeneous 
equation □5F^(X) = (with the Faddeev-Popov operator Ds). Therefore Fa{X) = everywhere in 
the bulk and the on-shcU contribution of the gauge-breaking term (4.4) vanishes and makes the whole 
on-shell action gauge independent. 
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over the fields on the brane to generate the full set of Feynman diagrams. This integral 
in which the bulk integration is done first, while the integration over brane fields is 
reserved for the last, schematically looks like 

Here eUipses denote all the details of gauge- fixing procedure (including the contribution 
of relevant Faddeev- Popov determinants), the subscript of D in the integration measure 
indicates the spacetime dimensionality of fields — path integral over the 5-dimensional 
metric Gab{X) in the bulk versus 4-dimensional fields g^v{x) and Na{x) on the brane. 
^\oop[g^lv^ Na] denotes the loop part of the brane effective action following from the 
functional integration over bulk fields bounded by brane fields g^u{x) and Na{x). 

In this form the functional integral is not manifestly covariant, because the inte- 
grand on the right-hand side is ofT shell and, therefore, gauge dependent. Its overall 
gauge independence is not manifest because it is attained only as a result of brane field 
integration. However, the strategy of its calculation can be improved by including the 
integration over D^^Na into the bulk, so that (6.6) takes the form 

j D^g^, e-^^n9,,^]-rxooM (...) (6.7) 

with the brane tree- level S^^[gfj_i,] and loop effective /loopb/ni/] actions gauge indepen- 
dent off shell. Such improvement, the corresponding Feynman diagrammatic technique 
and its gravitational Ward identities will be presented in the coming paper [41], and 
this technique is expected to establish in a simple way the covariant gauge-independent 
status of the strong-coupling problem*. 

To summarize, we extended the suggestion of [9] for infrared modifications of Ein- 
stein theory to the general case characterized by two nonlocal operators in Ricci tensor 
and scalar quadratic forms and considered restrictions on these operators imposed 
by the absence of VDVZ discontinuity. In models with low strong-couphng scale we 
showed that their effective action also contains the strongly coupled mode which cannot 
be perturbatively integrated out in terms of the metric variables. It serves in the lowest 
order as a Lagrange multiplyer for the constraint on metric variables which enforces the 
curvature scalar to vanish on shell. We demonstrated this situation on the examples 
of FP, DGP models and the recently suggested constrained version of the DGP model. 

^Depending on gauge-fixing procedure for both bulk and brane difi^eomorphisms the diagrammatic 
technique for this integral can be essentially simplified by factorizing the bulk and brane-to-brane 
parts [41] , which allows one efficiently to disentangle the different brane and bulk scales of the model. 
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A. Gravitational potentials in the constrained DGP 
model 

Bulk equations of motion for the action (4.1) modified by the term (5.1) read 

4'^ + ^^(M^^) = 0, (A.l) 



S^'J + -d,B, = 0, (A.2) 
^i? - -^d,B^ - 0, (A.3) 

where Sab = Rab — \GabR^^^ is the 5-dimensional Einstein tensor. Bianchi identities 
V^S^^ = then yield the equations of motion for gauge conditions 

= 0, (A.4) 
= -^d,^d,B,, (A.5) 

in which only is propagating in the bulk, while B^ express algebraically in terms of 
B5. 

Prom (A.2)-(A.3) one gets the 574-componcnts of the metric 



H,, = ^{d'^H,, - d,H) - (7 + 1) d,^^{ddH - UH), (A.6) 
1 1„1 



H,, = -ddH - ^P^iddH - DH), (A.7) 



whence 



^5,,,,, „ 1^, dl 



B, = d^^H,^ = -7 ^{ddH - OH), B^ = -pd^^iddH - □//), (A.8) 

where ddH = d^d'^Hj^^ and H = rj^^Hf^y. Combining (A.5), (A.8) and the definition 
of one has 

d^'B^ = D/iss - ddH = ^ S {^ddH -OH). (A.9) 
Together with (A.7) this implies 

n^iddH - an) = 0. (A.IO) 
The trace of (A.l) then gives 

n^H - (2-f + (3/2)(ddH -DH) ^0. (A.ll) 
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With the 4-dimcnsional transverse-traceless decomposition of H^i, {'q'^'^HjJ — 0, 



d'^HlJ = 0, d>^^, = 0) 



(A.12) 



the gauge becomes = —Oip^ + d^lh^^ — $ — D^). On the other hand, from (A. 5) 
and (A. 8) it equals i?^ = 3pd^ip/2. Comparison of these expressions gives </7^ = 0. 

Thus, finally from (A. 10), (A. 11) and the transverse-traceless part of (A.l) we have 
the set of bulk equations of motion for nonvanishing metric components 



and 



□si^JJ = 0, 
□5$ = 0, 

□5* + 3(27 + /5/2)$ = 



H,^ = 3^%^<!>, i/55 = (1 + 3/3/2)$ + 



(A.13) 
(A.14) 
(A.15) 

(A.16) 



Linearized Israel junction condition for the action (4.1) with the additional term 
(5.1) (in the presence of conserved matter sources on the brane) reads 



K^. - V,.K) - Ef} + ^(a^F, + d^F^ - v,.d^F") + |^T^, = 0, (A.17) 



MI 



where Kf^^ = {l/2)[d^H^y — d^H^.^ — d^H^^) is the extrinsic curvature of the brane. 
Transverse-traceless part of (A.17) gives 

IGtt, 



TT 



Ml '^^ ' 



while the longitudinal and trace parts imply 

1 



-□2l]> 



((97 - 6)m95 + - 3) °) * I - ^ (Srn^s + ^n) m 



Ml 



T = 0. 



(A.18) 

(A.19) 
(A.20) 



For /9 = 7 = the model essentially simplifies. Equations (A.14) and (A.15) 
decouple and under zero boundary conditions at y = 00 yield 



(A.21) 



so that everywhere in equations (A.18)-(A.20) ^5 can be replaced by the (nonlocal) 
operator on the brane —A. Their solution on the brane becomes 

Stt 1 „ ^ ^D-ScrmA 



3M| □ - 4mA - 3(7m2 



□2 



(A.22) 
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which finally leads to the gravitational potential (5.5)-(5.9). 

Note that the (redundant from the gauge-fixing viewpoint) brane gauge term in 
(A. 17) precludes from enforcing the bulk gauges and, thus, gives rise to gauge depen- 
dence discussed in Sect. 5. Taking the 4-dimensional divergence of the Israel junction 
condition (A. 17) and using the Gauss-Codazzi equation V^{K^y — ri^yK) = R^} one 
has in virtue of (A. 2) 



1 



7 



1 



-d^B^ = nF^. (A.23) 



ma 



Thus, the equation for gauge conditions in the bulk (A. 4) has nonzero data on the brane 
and has a nonvanishing solution. Therefore bulk gauge conditions are not enforced, 
and their contribution to equations of motion (A.1)-(A.3) (proportional to ^5/7) is 
nonvanishing even on shell. This situation holds also in the limit of 7 — > 0. 
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